Method and device for measuring the angular velocity of a luminance transition zone and steering aid system for fixation and tracking a target comprising at least one such luminance transition zone

ABSTRACT

The present invention relates to a method for measuring a time delay between two output signals measured by two adjacent optical sensors, the time delay resulting from a movement of an object moving in a given direction. The method comprises the steps:
         carrying out a spatial filtering of the signals Ph 1  and Ph 2  delivered by a first and a second optical sensor;   calculating a first order derivative and carrying out a temporal low-pass filtering of the signals;   calculating the second order derivative of the signals;   measuring the delay between the signals;
 
and is characterized in that the step of measuring the delay is a feedback loop based on an estimated delay between the temporally filtered signals.
       

     Device and set of devices for measuring the angular velocity of a luminance transition zone and steering aid system for fixation and tracking a target comprising at least one such luminance transition zone.

The present invention relates to a method and a device for measuring the angular velocity of a luminance transition zone of an object moving in a given direction and a system for fixing and tracking a target comprising at least one such luminance transition zone. For example, the object is a picture or a landscape having a luminance transition zone.

In particular, such steering aid systems for locating and tracking a target are used in microrobotics for the development of automatic flight control systems. For example, OCTAVE (for Optical Altitude Control system for Autonomous VEhicles) allowing a tethered aerial robot (the OCTAVE Robot) to perform terrain following and OSCAR (for Optical Scanning sensor for the Control of Autonomous Robots) allowing a tethered aerial robot (the OSCAR Robot) to locate and track a contrasting target with a high level of accuracy. OCTAVE and OSCAR use an optic flow sensor called EMD (for Elementary Motion Detector) that is described by:

F. Ruffier, S. Viollet, S. Amic and N. Franceschini in the article “Bio-inspired optic flow circuits for the visual guidance of micro-air vehicles”, IEEE Int. Symposium on Circuits and Systems ISCAS 2003, Bangkok, Thailand, pp. 846-849. This article is designated A₁.

F. Ruffier, S. Viollet, N. Franceschini in the article Visual control of two aerial micro-robots by insect-based autopilots Advanced Robotics, (2004) Vol. 18, No. 8, pp. 771-786. This article is designated A₂.

A₁ and A₂ describe a method and a device for measuring the angular velocity of a luminance transition zone that are based on an EMD (for Elementary Motion Detector) sensor comprising two optical sensors whose optical axes are angularly separated. The method described comprises the following steps:

-   -   a spatial filtering of the two optical sensors,     -   a band-pass temporal filtering of each signal delivered by the         two optical sensors,     -   thresholding each signal in order to measure the time difference         Δt between the signals delivered by the two optical sensors, the         time difference depending on the angular velocity of this         luminance transition zone.

However, the measurement of an angular velocity by the EMD sensory device can be disturbed by fast variations of luminosity, contrast, spatial frequency, linear velocity and distance from the objects in an unknown visual environment. Indeed, before thresholding, the amplitude of the visual signals varies according to the variations of these parameters. Because of the presence of low amplitude noise, it is difficult to set the value of the detection threshold below the chosen value. However, even if the detection threshold is set to an optimum level, the accuracy of the measurement of Δt will be influenced by the noise of the sensors present at that level, i.e. the noise present at the two intersections of the threshold, with the signals of the first and second optical sensors.

The output of the EMD is sometimes false, due to a mismatch, for example, when the two optical sensors detect at the same time not the same luminance transition zone, but two different ones. Due to the slow update rate of the EMD, the EMD holds a false measurement for a relative long time in this case. The slow update rate of the EMD generates a discontinuous signal composed of transient signals (stairs) that may cause many problems when used in control loops.

Moreover, the output of the EMD is not directional and gives even wrong results when the EMD sensor measures angular velocities in the direction opposed to its preferred direction.

The aim of the invention is to solve this problem by modifying the signal processing in order to obtain a more robust, bidirectional measurement of the angular velocity of luminance transition zones of the visual scene with a high refresh rate.

Therefore, the invention relates to a method for measuring a time delay between two output signals measured by two adjacent optical sensors, the time delay resulting from a movement of an object moving in a given direction, the method comprising:

carrying out a spatial filtering of the signals Ph₁ and Ph₂ delivered by a first and a second optical sensor, the optical axes of the first and the second optical sensors forming an inter-receptor angle (Δφ) containing the object;

calculating a first order derivative and carrying out a temporal low-pass filtering of the signals delivered by the first and second optical sensors in order to obtain signals Ph₁′ (t) and Ph₂′ (t);

calculating the second order derivative of the signals delivered by the first and second sensors in order to obtain signals Ph₁″ (t) and Ph₂″ (t);

measuring the delay (Δt) between the signals delivered by the first and second optical sensors, the delay depending on the angular velocity of the luminance transition zone of the object;

characterized in that the step of measuring the delay Δt is a feedback loop based on an estimated delay between the temporally filtered signals delivered by the first and second optical sensors, the feedback loop comprising the steps of:

estimating a time delay Δt(t_(n)) at an instant t_(n) between the temporally filtered signals delivered by the first and second optical sensors;

calculating the error Ph_(diff)(t_(n)) at the instant t_(n) between the temporally filtered signal delivered by the second optical sensor Ph₂′ (t_(n)) at the instant t_(n) and the temporally filtered signal Ph₁′ (t_(n)−Δt(t_(n))) delivered by the first optical sensor and delayed by the estimated delay Δt(t);

estimating whether Ph_(diff)(t_(n))≠0 i.e. whether a delay error Δt_(error)(t_(n)) exists;

if Ph_(diff)(t_(n))≠0, estimating whether the delay error Δt_(error)(t_(n)) is smaller or bigger than zero on the basis of the error Ph_(diff)(t_(n)) and the second order derivative Ph₁″(t) and Ph₂″(t) of the signals delivered by the first and second sensors; and

calculating the evolution of the delay Δt(t) at the instant t_(n) on the basis of Δt(t_(n)) and Δt_(error)(t_(n)).

In other embodiments, the method comprises one or several of the following features, taken in isolation or in any technically feasible combination:

the feedback loop on the estimated delay comprises the steps of:

-   -   filtering the temporally filtered signals of the first optical         sensor and/or the second optical sensor with a derivative filter         in order to compute the second-order temporal derivative         function Ph₁″(t) and/or Ph₂″(t) of the first optical sensor         and/or second optical sensor at an instant t_(n);     -   determining the temporal evolution of the estimated delay Δt_(n)         on the basis of the error Ph_(diff)(t_(n)) and the second-order         temporal derivative function Ph₁″(t) or Ph₂″(t) and updating the         estimated delay Δt     -   whereby these steps are continuously repeated until the         estimated delay error Δt_(error)(t) between the delayed signal         of the first optical sensor and the signal of the second optical         sensor is substantially zero;

the update of the estimated delay is increased or decreased with a constant slope according to the sign of the estimated delay error (Δt_(error)(t_(n))) which is the ratio of the error Ph_(diff)(t_(n)) to the second-order temporal derivative function Ph₁″(t_(n)) or Ph₂″(t_(n));

the update of the estimated delay is increased or decreased by the product of a function α(t);

the error Ph_(diff)(t) between the delayed signal of the first optical sensor and the signal of the second optical sensor or between the delayed signal of the second optical sensor and the signal of the first optical sensor reflects the accuracy of the time delay estimation;

the function α(t) is a constant value between 0 and 1;

the function α(t) is proportional to the square of the estimated delay Δt;

the function α(t) is modulated by the function |cos(Δt*2πF)|, where Δt is the estimated delay and F is a frequency present on background signals;

the feedback loop on the estimated delay comprises the step of resetting and initializing an integrator by a threshold method applied to the temporally filtered signals delivered by the first and second optical sensors if the error Ph_(diff)(t) exceeds a predetermined value;

the method comprises a step of detecting the direction of the movement of the luminance transition zone comprising the steps of

-   -   calculating a) Ph_(diff)(t) according to         Ph_(diff)(t)(A)=Ph₂′(t−Δt)−Ph₁′(t) and b) Ph_(diff)(t) according         to Ph_(diff)(t)(B)=Ph₁′(t−Δt)−Ph₂′(t)     -   comparing the temporal evolution of Ph_(diff)(t)(A) obtained         according to a) to Ph_(diff)(t)(B) obtained according to b)     -   determining on this basis whether the luminance transition zone         moves from the first to the second optical sensor or in the         opposite direction.

The invention also relates to a device for measuring an angular velocity of a luminance transition zone of an object moving in a given direction, the device comprising at least:

a first and a second optical sensor having each an optical axis (O₁Y₁, O₂Y₂) that delimit an inter-receptor angle (Δφ), Δφ containing the object;

a lens in order to carry out a spatial filtering of the signals delivered by the first and second optical sensors, the lens having an optical center, the first and second optical sensors being placed substantially in the image focal plane of this lens, the two lines linking the centers (O₁, O₂) of the first and second optical sensors and the optical center of the lens being defined as average directions of observation of the first and second optical sensors;

a detection circuit connected to the first and second optical sensors and comprising a band-pass filter in order to carry out a band-pass temporal filtering of the signals delivered by the first and second optical sensors, a means for measuring the delay Δt between the signals delivered by the first and the second optical sensors, the delay depending on the angular velocity of the light transition zone, and a calculator for calculating the angular velocity starting from the delay Δt, the angular velocity being the ratio of the inter-receptor angle (Δφ) to the delay Δt;

characterized in that means for measuring the delay Δt comprises:

a means to sample the temporally filtered signals of the first and second optical sensors, this means following the band-pass filter;

a means to estimate a delay Δt(t_(n)) at an instant t_(n) between the temporally filtered signals delivered by the first and second optical sensors;

a means to compute the error Ph_(diff)(t_(n)) at the instant t_(n) between the temporally filtered signal Ph′₂(t_(n)) delivered by the second optical sensor at the instant t_(n) and the temporally filtered signal delivered by the first optical sensor and delayed by the estimated time delay Ph₁′(t_(n)−Δt(t_(n))) at the instant t_(n); and

a means to compute the evolution of the estimated delay Δt(t) starting from the estimated delay Δt(t_(n)) at the instant t_(n) and the error Ph_(diff)(t_(n)) calculated at the instant t_(n).

In other embodiments, the device for measuring an angular velocity comprises one or several of the following features, take in isolation or in any technically feasible combination:

the first and second optical sensors are each formed by a photoelectric sensor;

the means for measuring the delay Δt comprises a derivative filter in order to calculate the second-order temporal derivative function Ph₁″(t) and/or Ph₂″(t) of the first optical sensor and/or of the second optical sensor at an instant t_(n) and a means to determine the temporal evolution of the delay Δt on the basis of the error Ph_(diff)(t) and the second-order temporal derivative function Ph₁″(t) and/or Ph₂″(t) and in order to update the estimated delay;

the detection circuit is a bidirectional detection circuit adapted to detect the direction of the motion of the luminance transition zone, that is, from the first to the second optical sensor or from the second to the first optical sensor, the bidirectional detection circuit comprising two means for measuring the delay Δt between the signals.

The invention also relates to a set of devices for measuring an angular velocity comprising at least one lens, a plurality of at least three optical sensors and a plurality of detection circuits, characterized in that two adjacent optical sensors form a pair and are placed in the focal plane of one of the at least one lens and connected to one of the plurality of detection circuits; in that each lens of the at least one lens and all optical sensors are arranged substantially on a spherical or cylindrical surface surrounding a common centre, the common centre constituting, the optical centre of an inter-receptor angle which is the sum over all inter-receptor angles of all elementary devices taken together, and in that the devices are devices according to the invention.

In other embodiments, the set of devices for measuring an angular velocity comprises one or several of the following features, taken in isolation or in any technically feasible combination:

the at least three optical sensors are arranged in two directions substantially orthogonal;

the optical sensors have a spectral sensitivity which can be the same or different for the individual optical sensors.

The method and device for measuring the angular velocity of a luminance transition zone based on a time delay measurement, according to the present invention, are used for industrial implementation of automatic steering aid systems and collision avoidance system.

The invention also relates to a steering aid system for the visual fixation and fine tracking of a target, the target being an object comprising at least one contrasting edge having a luminance transition zone, and for controlling the speed of an aircraft, characterized in that it comprises:

-   -   a device for detecting the angular position of the contrasting         edge;     -   a means to maintain the line of sight of the steering aid system         constantly on the position of the contrast edge;     -   a device or a set of devices for measuring an angular velocity         of the luminance transition zone according to the invention; and     -   a means to control the speed of the aircraft according to the         angular position and angular velocity of the target.

The invention will be better understood from the description that follows, which is provided solely by way of example and which refers to the drawings, wherein:

FIG. 1 is a schematic representation of the device for measuring an angular velocity of a luminance transition zone according to the invention;

FIG. 2 is a schematic representation of a part of the detection circuit of the device for measuring an angular velocity according to the invention;

FIG. 3 is a schematic representation of the method for measuring the angular velocity of a luminance transition zone according to the invention;

FIG. 4 is a schematic representation of the directional sensitivity of the optical sensors of the device for measuring the angular velocity according to the invention;

FIG. 5 is a schematic representation of an example of the implementation of the method for measuring the angular velocity of a luminance transition zone according to the invention;

FIGS. 6 to 9 are schematic representations of different embodiments of an estimator of the time delay between the two signals delivered by two optical sensors of the device for measuring an angular velocity according to the invention;

FIG. 10 is a schematic representation of another embodiment of the device for measuring an angular velocity which is bidirectional according to the invention;

FIG. 11 is a schematic representation of the implementation of the method for measuring the angular position of a contrasting edge according to the invention;

FIG. 12 is a schematic representation of the overall device for measuring the angular position of a contrasting edge according to the invention;

FIG. 13 is a schematic representation of the directional sensitivity of the optical sensors integrated in the device for measuring the angular position of a contrasting edge according to the invention;

FIGS. 14 and 15 represent the calculated device output signal (Y_(out)(t)) as a function of the angular position Ψ_(C) of a contrasting edge (FIG. 14) and of a contrasting bar (FIG. 15) according to the invention;

FIG. 16 represents the calculated device output signal (Y_(out)(t)) as a function of the angular position Ψ_(C) of a contrasting bar, plotted for various ratios of the full width at half maximum (FWHM) Δρ to the inter-receptor angle Δφ;

FIG. 17 shows the signals obtained at the different steps of the method for measuring the angular position according to the invention, when considering a sinusoidal vibrational law whose amplitude (0.25°) is much smaller than the inter-receptor angle Δφ,

FIG. 18 is a schematic representation of another embodiment of the device for measuring the angular position of a contrasting edge according to the invention;

FIGS. 19 and 20 are schematic representations of yet another implementation of the method for measuring the angular position of a contrasting edge according to the invention; and

FIGS. 21 to 23 are schematic representations of different embodiments of a demodulator integrated into the device for measuring the angular position of a contrasting edge according to the invention.

A more detailed description of the device and the method for measuring the angular velocity of a luminance transition zone which is substantially rectilinear in a given direction, according to the present invention, will now be given with reference to the FIGS. 1 to 10.

Referring to FIG. 1, the luminance transition zone and, in particular, a contrasting edge E of an object P are rectilinear in a given direction, along the Z axis which is substantially orthogonal to the plane containing FIG. 1.

The device 2 according to the invention comprises at least one optical device 4 comprising a first and a second optical sensor D₁ and D₂, separated by an inter-receptor angle, denoted by Δφ. The inter-receptor angle Δφ is the angle between the two optical axes O₁Y₁ and O₂Y₂ of the first and second optical sensors D₁ and D₂, respectively, with O₁ and O₂ being the centers of the first and second optical sensors, respectively.

The object P and its luminance transition zone are moving with a velocity denoted by V which is represented by an arrow in FIG. 1. The Object P is moving along the X axis which is orthogonal to the given direction along the Z axis. For example, object P is moving in the direction from the first to the second optical sensor. Consequently, luminance transition zone E is first detected by the first optical sensor D₁ and then by the second one D₂.

In a non-limiting embodiment, the optical device 4 comprises a convex lens L which determines this inter-receptor angle Δφ. The optical sensors D₁ and D₂ are placed substantially in the image focal plane of the lens L, so that the optical center O of the lens L is situated between the optical sensors D₁, D₂ and the luminance transition zone.

A reference direction is advantageously the direction OY₁₂ substantially corresponding to the average direction of the angle bisector between the optical axes O₁Y₁ and O₂Y₂ of the first and second optical sensors D₁, D₂.

In a preferred non-limiting embodiment, the first and the second optical sensors D₁ and D₂ are formed by photoelectric sensors such as photoelectric diodes, whose spectral sensitivity can be selected either in the spectral range of visible or even ultraviolet light or, alternatively, in the range of near of far infrared radiation, for example for nocturnal detection.

The wavelength of maximum sensitivity of the aforementioned photodiodes can thus be chosen in dependence of the application of the device for detecting a contrasting edge E according to the present invention.

Moreover, the device 2 for measuring the angular velocity of the contrasting edge E according to the invention comprises a detection circuit 6 in order to compute, on the basis of the signals delivered by the first and the second optical sensors D₁ and D₂, the angular velocity ω(t) of the luminance transition zone.

The detection circuit 6 is connected to the first and second optical sensors and receives the signals delivered by the first and second optical sensors D₁ and D₂ which are denoted by Ph₁(t) and Ph₂(t), respectively.

The detection circuit 6 comprises a band-pass filter 8 connected to the optical sensors in order to compute the temporal derivative signals of Ph₁ and Ph₂, denoted by Ph₁′ and Ph₂′, of the first and second optical sensors, respectively.

Furthermore, the detection circuit 6 comprises a means 10 for measuring the time delay Δt between the signals delivered by the first and the second optical sensors, the delay depending on the angular velocity of the luminance transition zone.

Optionally, the detection circuit 6 can comprise a calculator 12 for calculating, for example the angular velocity ω(t) (optic flow), or any mathematical function of the delay Δt (calculated between the output signals of the optical sensors). The angular velocity ω(t) is the ratio of the inter-receptor angle Δφ to the delay Δt:

${\omega (t)} = \frac{\Delta\phi}{\Delta \; {t(t)}}$

Means 10 for measuring the delay Δt will now be described in detail with reference to FIG. 2. Means 10 forms a feedback loop 13 in order to determine the temporal evolution of the delay Δt and therefore to update the delay with every cycle of the feedback loop.

The feedback loop 13 for measuring the delay Δt comprises at least means 14 to delay the temporally filtered signal Ph₁′ delivered by the first optical sensor D₁, an estimator 16 of the temporal error Δt_(error) of the delay Δt and a means 18 to set the delay Δt to a value Δt_(init) calculated by an EMD when an error Ph_(diff), which is initializing means 18, is above a predetermined limit during a predetermined time.

The means 10 for measuring the delay Δt is shown in FIGS. 2 and 6. In this embodiment means 10 comprises:

-   -   means 13 to update the estimated delay Δt(t) at an instant t_(n)         between the temporally filtered signals delivered by the sampled         first and second optical sensor output signals, starting from         the estimated delay Δt(t_(n)) and the temporal error         Δ_(error)(t) at the instant t_(n);     -   a means 15 to compute an error Ph_(diff)(t) at an instant t_(n)         between the temporally filtered signal delivered by the second         sensor at the instant t_(n) Ph₂′(t_(n)) and the temporally         filtered signal delivered by the first optical sensor at the         instant t_(n) delayed by the estimated delay Δt(t) at the         instant t_(n):Ph₁′(t_(n)−Δt(t_(n)));     -   a derivative filter 17 in order to compute the second-order         temporal derivative function of the second optical sensor at the         instant t_(n):Ph₂″(t_(n)); and     -   a means 19 to determine whether the estimated delay Δt(t) at the         instant t_(n) should be increased or decreased in order to         approach the real value of Δt(t).

A more detailed description of the method 100 for measuring the angular velocity of the luminance transition zone according to the present invention will now be given with reference to FIGS. 1 and 3 to 5.

When a luminance transition zone moves in front of the optical device 4, along the X axis in the direction from the first optical sensor D₁ to the second one D₂, the sensors D₁ and D₂ convert the luminance coming from the object, and in particular from the luminance transition zone, in analogical signals (current or voltage). The temporal signals delivered by the first and second sensors D₁ and D₂ are denoted by Ph₁(t) and Ph₂(t), respectively.

The method 100 according to the invention comprises a step 102 of spatial filtering the signals delivered by the optical sensors D₁ and D₂.

This spatial filtering 102 is determined by the lens L (FIG. 3). Indeed, the lens confers to the two optical sensors a Gaussian sensitivity function, i.e. a Gaussian angular directivity, as shown in FIG. 4. These two directional sensitivity functions of the two optical sensors D₁ and D₂ overlap partially. The Gaussian shape is due to the convex lens defocusing the incoming light. The maximum sensitivity of each of the optical sensors D₁ and D₂ is along the optical axes O₁Y₁ and O₂Y₂ which delimit the inter-receptor angle Δφ.

The first and second optical sensors D₁ and D₂ have a directional sensitivity, according to the angular position Ψ of the light source:

${s(\psi)} = {\frac{1}{\sigma \sqrt{2\pi}}^{- \frac{\psi^{2}}{2\sigma^{2}}}}$

The parameter Δρ shown in FIG. 4 is the full width at half maximum (FWHM) of the sensitivity function and results directly from the variance σ² of the Gaussian function with Δρ=σ×2√{square root over (2 ln 2)}.

Typical optical characteristics of the optical device 4 are: Δρ=3° and Δφ=2.87°.

Placed in front of a contrasting edge, each optical sensor measures a value that is the integral of the angle ψ, of the sensitivity function s(ψ) of the optical sensor multiplied by the luminance of the contrasting edge E.

An example of the spatially filtered signals is shown in FIG. 5( a). The signals Ph₁(t) and Ph₂(t) have the same shape but Ph₂(t) is time delayed by Δt compared to Ph₁(t) due to the inter-receptor angle Δφ between the two optical axes of the optical sensors.

Then, in step 104 (FIG. 3) a band-pass filter 8, which is generally an analogical circuit, temporally filters the signals Ph₁(t) and Ph₂(t) delivered by the first and second optical sensors D₁, D₂ in order to compute the temporally derivative signals of the optical sensors, which are denoted by Ph₁′(t) and Ph₂′(t) (FIG. 5 (b)). The band-pass filter 8 comprises a high-pass filter and a low-pass filter.

The high-pass filter acts as a temporal differentiator of the signals delivered by the first and second optical sensors D₁ and D₂, especially in the frequency range from 0 to 30 Hz which contains the relevant information. The high-pass filter also cancels the DC component of the signals Ph₁(t) and Ph₂(t).

The low-pass filter attenuates the high frequency components of the luminance and reduces noise and interferences.

After temporal filtering 104, the temporally differentiated signals are sampled, for example at a sampling frequency of 2 kHz.

In step 106, means 10 determines the delay Δt between the signals Ph₁(t) and Ph₂(t) and in step 108, the calculator 12 calculates the angular velocity ω(t) according to the equation:

${\omega (t)} = \frac{\Delta\phi}{\Delta \; {t(t)}}$

or another function of the delay Δt.

The determination 106 of the delay Δt will now be described in detail. The delay Δt(t) is determined in a feedback loop The feedback loop comprises the steps of:

(a) estimating the delay Δt(t) between the temporally filtered signals Ph₁′(t) and Ph₂′(t) at an instant t_(n), this delay being denoted by Δt(t_(n));

(b) delaying the temporally filtered signal Ph₁′(t) of the first optical sensor by the delay Δt(t_(n)). The delayed signal is denoted by Ph₁′(t−Δt(t_(n))). An example of this step is represented in FIGS. 5( c) and 5(d), which is the enlarged view of the part shown in a box in FIG. 5( c);

(c) calculating the error Ph_(diff)(t) at the instant t_(n) between the temporally filtered signal Ph₂′(t) delivered by the second optical sensor at the instant t_(n) and the temporally filtered signal Ph₁′(t−Δt(t_(n))) delivered by the first optical sensor and delayed by the estimated delay Δt(t_(n)) at the instant t_(n):

Ph _(diff)(t _(n))=Ph ₁′(t _(n) −Δt(t))−Ph ₂′(t _(n))

If the estimated delay Δt(t_(n)) is correct, i.e. corresponds to the real delay Δt at the instant t_(n), then Ph_(diff)(t_(n)) is zero. Otherwise Ph_(diff)(t_(n)) is different from zero.

When the output of the feedback loop Δt(t_(n)) is incorrect, the feedback loop takes a signal Δt_(error)(t) as an error signal which will be integrated and added to the current value of Δt(t) in order to correct Δt so that Δt converges toward the real value of Δt(t) and that the error Ph_(diff)(t) converges toward zero.

The determination of the error Δt_(error)(t) shall now be considered in detail.

In FIG. 5( d), the signals Ph₂′(t) and Ph₁′(t−Δt) have a negative slope at the instant t_(n) (Ph₂″(t_(n))<0 and Ph₁″(t_(n)−Δt_(n))<0). The error Ph_(diff)(t_(n)) according to the equation above is also negative. As it can be seen from FIG. 5( d), the estimated delay Δt is too small. Consequently, Δt has to be augmented by a positive value Δt_(error)(t_(n)) in order to reduce Ph_(diff)(t). It is important to note that in the case (not shown) that the signals Ph₂′(t) and Ph₁′(t−Δt) have a positive slope at the instant t_(n), and that Ph₁′(t−Δt) is larger than Ph₂′(t), the error Ph_(diff)(t_(n)) would be positive and the estimated delay Δt too small. Consequently, the value Δt_(error)(t_(n)) needs to be positive in order to reduce Ph_(diff)(t).

In the opposite case (not shown), that is when the estimated delay Δt(t_(n)) is too large, Ph_(diff)(t_(n)) is smaller than zero for signals with a positive slope at the instant t_(n) and larger than zero for signals with a negative slope at the instant t_(n). Consequently, the direction of evolution of the estimated delays Δt(t) at the instant t_(n) depends on the one hand on the sign of the slope of Ph₂′(t) (or of Ph₁′(t−Δt) at the instant t_(n) (sign of Ph₂′(t) or of Ph₁′(t_(n)−Δt)) and on the other hand on the error Ph_(diff)(t) at the instant t_(n).

In the case that the error Ph_(diff)(t_(n)) is different from zero, the evolution of the delay Δt(t) at the instant t_(n) should comply with the correlation given in the following table:

The feedback loop comprises furthermore the steps of:

(d) filtering the filtered signal of the second optical sensor Ph₂′(t) by a derivate filter in order to compute the second-order temporal derivative function Ph₂″(t) of the signal of the second optical sensor at the instant t_(n), denoted by Ph₂″(t_(n));

(e) estimating the temporal error Δt_(error)(t_(n)) of the delay Δt(t_(n)), on the basis of the error Ph_(diff)(t_(n)) and the second-order temporal derivative function Ph₂″(t_(n)) (or the signal Ph₁″(t_(n)−Δt(t_(n)))); and

(f) estimating the evolution of Δt(t) from the temporal error Δt_(error)(t_(n)). Many alternatives exist. A first one consists in setting a fixed slope value, denoted by α, of the evolution, and then using the sign of Δt_(error)(t_(n)) to define the direction of evolution at each instant. The value of Δt(t) is then updated according to:

Δt(t)=Δt _(init)+∫_(init) ^(n)α×2*(sign(Δt _(error)(t))−0.5)dt

wherein the sign function takes the value+1 for positive values of Δt_(error)(t), and 0 for negative values. Other alternatives consist in defining a slope α proportional to the absolute value Δt_(error). This means that α becomes a function of the error Δt_(error)(t).

Moreover, the feedback loop can reinitialize the value of Δt to a value Δt_(init). This reinitialization occurs when the error Ph_(diff)(t) between the delayed signal of the first optical sensor and the signal of the second optical sensor is important. The reinitialization of Δt may be triggered when the error Ph_(diff)(t) stays above a predetermined value during a predefined time.

The reinitialization value Δt_(init) is the output of a classical Elementary Motion Detector (EMD) as described in article A₁.

Now, four alternative methods, represented in FIGS. 6 to 9, will be described. They use different estimators 16 of the delay.

FIG. 6 shows the first embodiment of the estimator 16. The estimated delay Δt(t) increases or decreases with respect to the delay Δt(t_(n)) with constant slope in dependence of the sign of the product of the error Ph_(diff)(t_(n)) and the second-order temporal derivative function Ph₂″(t_(n)).

If the error Ph_(diff)(t_(n)) has the same sign as the second-order temporal derivative function Ph₂″(t_(n)), Δt(t_(n)) is too small and Δt(t) has to increase with slope ε. If Ph_(diff)(t_(n)) and Ph₂″(t_(n)) have opposite signs, Δt(t_(n)) is too large and Δt(t) has to decrease with slope ε.

A simplified expression determining the evolution of Δt(t) is

Δt(t)=Δt _(init)+∫_(init) ^(n)×2*(sign(Ph _(diff)(t))⊕sign(Ph″ ₂(t))−0.5)dt

where ε is a predetermined value and ⊕ is the exclusive or operator (xor). The sign function takes the value+1 for positive values, and 0 for negative values. This method is extremely simple to implement as the functions ⊕ (xor) and sign are basic and belong to the functions that are most quickly carried out in a microcontroller.

However the speed of evolution of the delay depends on the slope of +ε or −ε.

A second alternative of the estimator 16 shown in FIG. 7, takes into account the values of the error Ph_(diff)(t) and of the second-order temporal derivative function Ph₂″(t). The temporal error can be expressed according to the following approximation (FIG. 5( d)):

${\Delta \; {t_{error}\left( t_{n} \right)}} \approx \frac{{Ph}_{diff}\left( t_{n} \right)}{{Ph}_{2}^{''}\left( t_{n} \right)}$

In practice, the signals are noisy, in particular Ph₂″(t). Also, this relation for the evolution for Δt used as such gives a disturbed result and in the worst case it cannot lead to the solution. So, the ratio is weighted/ by a factor α chosen in the interval [0; 1]:

${\Delta \; {t_{error}\left( t_{n} \right)}} = {\alpha \times \frac{{Ph}_{diff}\left( t_{n} \right)}{{Ph}_{2}^{''}\left( t_{n} \right)}}$

The smaller the factor α, the smoother will be the convergence of the estimated delay Δt. However, a small value of α will decrease the speed of convergence of the estimated delay Δt, which complicates the estimation of the delay, provided that Δt varies quickly (for example, during a rapid change of the angular speed rate).

In the two preceding alternative methods, the measurement noise on the estimation of Δt is constant and does not depend on the value of Δt itself. Thus, the computation of the angular velocity, which is reciprocally proportional to Δt, makes this noise not constant on the whole measurement range of ω(t). The measurement noise will be higher for high angular speed. The same reasoning can be applied for the speed of convergence of the method for computing Δt. This convergence speed or dynamic is constant over the whole range of the Δt measurement. Thus, the dynamic on the resulting angular speed ω(t) depends on the angular speed itself i.e. the higher the angular speed, the faster the dynamic.

A third embodiment of the estimator 16 shown in FIG. 8, takes into account the values of Δt to estimate the error Δt_(error). The corresponding method uses, instead of a factor, a function of the angular velocity ω(t) (or of its inverse Δt)). This function, denoted by α(Δt(t), is chosen such as to get a constant noise on the computed angular velocity ω(t). The estimated value of an angular velocity variation will then have the same dynamic within the entire measurable range of angular velocity.

The update of the estimated delay is increased or decreased by the product of a function α(t) and the ratio of the error Ph_(diff)(t) to the second-order temporal derivative function Ph₂″(t). The function α(Δt(t) is proportional to the square of the estimated delay according to the equation:

α_(a)(Δt(t))=α₀×(Δt(t))²

wherein α₀ is a factor chosen in the interval [0; 1]. In a fourth embodiment, the function α(Δt(t_(n))) is modulated by |cos(2πΔt/10e⁻³)|:

α_(a)(Δt(t _(n)))=α₀×(Δt(t _(n)))²×|cos(2πΔt/10e ⁻³)|

The corresponding estimator 16 is shown in FIG. 9. The temporal error can be calculated according to the following equation:

α_(a)(Δt(t))=α₀×(Δt(t))²×cos(2πΔt(t)/10e ⁻³)

When the sensor is used in the presence of an artificial light (neon, 100 Hz), the artificial light is mainly filtered by the analogical filters of the step (a). However, when weak angular velocities are measured, the amplitude of the temporally derived function of the signals is low and the residue at 100 Hz, which does not contain relevant information, becomes a possible source of noise on the sensor output. This residue at 100 Hz has a weak influence when the delay Δt is a multiple of 10 ms because the comparison of Ph₂′ and delayed Ph₁′ is carried out in phase with this noise at 100 Hz. Therefore, in this alternative method, α(Δt) has to be modulated by the function |cos(Δt*2πF)| where F is the 100 Hz frequency.

In another embodiment of the device 2 for measuring the angular velocity according to the invention and shown in FIG. 10, the device 2 comprises a bidirectional detection circuit 200 which comprises a means 202 to detect the direction of the movement of the luminance transition zone E.

The bidirectional detection circuit 200 comprises, as the aforementioned detection circuit 6 (FIG. 1), at least a band-pass filter 8 connected to the optical sensors in order to compute the derivative signals, denoted by Ph₁′ and Ph₂′, of the first and second optical sensors, and a calculator 12 (not shown) for calculating the angular velocity ω(t).

Furthermore, the bidirectional detection circuit 200 comprises a first and second means 10-A and 10-B for measuring the delay Δt between the signals delivered by the first and second optical sensors, each of these means corresponding to means 10 (FIG. 1) for measuring the delay Δt according to the invention and as described before. Each of means 10-A and 10-B for measuring the delay is connected between the band-pass filter 8 and the calculator 12.

Means 10-B delays the signal delivered by the first optical sensor Ph₁′(step (b)) and uses the signal Ph₂′ as such. In contrast to means 10B, means 10-A delays the signal delivered by the second optical sensor Ph₂′ and uses the signal Ph₁′ as such.

In the following, the same steps (c) to (f) as described before are executed. Means 202 which detects the direction of the movement of the luminance transition zone is connected between means 10-A and 10-B on the one hand, and the calculator 12 on the other hand.

Then, means 202 compares the errors Ph_(diff)(A) and Ph_(diff)(B) calculated by the first and second means 10-A and 10-B, respectively, and uses the output of that means 10-A or 10-B for which Ph_(diff) converges to zero.

Indeed, means 202 selects first means 10-B. In that case, means 202 can deduce the angular velocity if the luminance transition zone moves from the first to the second optical sensor. And then means 202 selects means 10-A. In this case, means 202 can deduce if the luminance transition zone moves from the second to the first optical sensor.

A device 2 comprising such a bidirectional detection circuit 200 is a bidirectional device because it works in the direction of the displacement of the luminance transition zone of the object, as well as in the opposite direction.

Measurements have shown that the dynamics of the detection device according to the invention, i.e. the capacity to measure an abrupt variation of the measured value, is as fast as the “classical” device.

The detection device for measuring an angular velocity according to the invention brings the following advantages. According to one embodiment of the invention, the measurement of the angular velocity can be bidirectional. Experiments have shown that this device can be used in a wide range of luminosity. The measurement is not disturbed by a variation of luminosity or a sudden variation of the speed of the optical device relative to its visual environment (absence of transitory time). The angular velocity data are continuously updated. The device is simple and requires only few resources of calculation of a microcontroller. Its implementation is fast and measurements with the method are robust and more precise than by a classical method using thresholding.

According to an embodiment, several devices for measuring the angular velocity according to the invention can be combined forming a set of devices. Such a set of devices comprises at least one lens, a plurality of N optical sensors (N being an integer and N≧2) placed in the focal plane of the lens and a plurality of N−1 detection circuits as described hereinbefore. This embodiment has an increased total field of vision, which can even be panoramic, if necessary. Each pair of adjacent optical sensors is connected to a detection circuit to compute a plurality of output signals representing the angular velocity ω(t) within each inter-receptor angle of two adjacent optical sensors. The output signals delivered by the N−1 detection circuits can thus be sampled for subsequent digital processing. The method as described hereinbefore is applied to each pair of optical sensors to measure the angular velocity of a luminance transition zone.

The plurality of N optical sensors can be placed in the focal plane of a single lens, or each pair of consecutive sensors can be placed in the focal plane of one lens among a plurality of lenses.

Each lens and each optical sensor are arranged substantially on a spherical or cylindrical surface surrounding a common centre. This common centre constitutes, for the set of devices, the optical centre of an inter-receptor angle Δφ being the sum over all inter-receptor angles Δφ_(i) of each of the devices taken together.

Furthermore, this set of devices can be generalized for detecting a two-dimensional target comprising at least two contrasting edges (E), substantially rectilinear and orthogonal to each other.

In the following example, the number of optical sensors is limited to three, D₁, D₂ and D₃, in order to facilitate the description. The two pairs of optical sensors, D₁-D₂ and D₂-D₃, are arranged in two directions H and V substantially orthogonal to each other. Each pair of optical sensors is placed in the focal plane of one of two lenses or all the optical sensors are placed in the focal plane of a single lens. Each lens and each optical sensor are arranged substantially on a spherical or cylindrical surface surrounding at least one common centre. A first centre constitutes the optical centre of an inter-receptor angle Δφ_(H), in the first direction, and a second centre, which can be identical to the first one, constitutes the optical centre of an inter-receptor angle Δφ_(V), in the second direction, substantially orthogonal to the first one.

In the same manner as before, each pair of optical sensors is connected to a detection circuit to compute the angular velocities ω_(H) (t) and ω_(V) (t). Each output signal depends on the position of the contrasting edge (E) of the target.

Finally, in a non-limiting embodiment of the set of devices, the plurality of optical sensors is arranged according to a matrix of optical sensors, each optical sensor constituting a pixel D_(ij). The rectangular matrix of optical sensors made up in this way forms a two-dimensional retina.

A more detailed description of a particularly advantageous embodiment of the device for measuring the angular velocity according to the present invention will now be given. In general, the reference direction OY₁₂ of the device can be any direction. In a non-limiting embodiment, however, this direction can advantageously be defined as perpendicular to the direction in which the luminance transition zone, and so the contrasting edge E, moves and to the direction in which it extends.

The corresponding method may be facilitated by orienting the reference direction OY₁₂ orthogonally to a plane substantially corresponding to the focal plane of the lens L.

For this purpose and in accordance with an embodiment of the detection device according to the present invention, the detection device advantageously comprises a means to orient the optical device 4 and with it the reference direction OY₁₂ of the detection device in the direction transverse to said direction in which the substantially rectilinear luminance transition zone extends and to said direction in which the luminance transition zone moves.

Such means for orientation are described in patent document FR 2 869 494.

A more detailed description of a steering aid system for the visual fixing and fine tracking of a target comprising at least one contrasting edge (E) having a luminance transition zone and for controlling the speed of an aircraft will now be given.

The steering aid system aims at the visual fixing and fine tracking of a target and at the control of the speed of an aircraft. The target comprises at least one contrasting edge having a light transition zone.

The system comprises a device for measuring the angular velocity of the luminance transition zone according to the invention, a device for detecting the angular position of the contrasting edge (E), a means to maintain the line of sight of the steering aid system constantly on the position of the contrasting edge E and a means to control the speed of the aircraft according to the angular position of the target and its angular velocity relative to the speed of the aircraft.

In the following, a detailed description of the device for detecting the angular position of the contrasting edge E and of a corresponding method will be given with reference to FIGS. 11 to 23.

Referring to FIG. 11, the luminance transition zone and, in particular, the contrasting edge E of an object P are rectilinear in a direction along the Z axis which is substantially orthogonal to the plane containing FIG. 11.

In reference to FIG. 12, the device for measuring the angular position according to the invention comprises at least a first and a second optical sensor D₁ and D₂ whose optical axes are separated by an inter-receptor angle Δφ.

In a non-limiting embodiment, this inter-receptor angle Δφ is determined by a convex lens L, as shown in FIG. 11, the optical sensors D₁ and D₂ being placed substantially in the focal plane of the lens L. The optical center O of the lens L is located between the optical sensors D₁, D₂ and the luminance transition zone. The optical axes of the first and second optical sensors D₁, D₂ correspond substantially to the lines linking the centers O₁ and O₂ of the first and second optical sensors, respectively, and the nodal point of the lens assembly, which for the sake of simplicity shall be called here the optical center O of the lens L. The optical axes O₁Y₁ and O₂Y₂ delimit the inter-receptor angle denoted by Δφ=(O₁Y₁, O₂Y₂)=(OO₁, OO₂).

The angular position of the contrasting edge E relative to a reference direction lying within the inter-receptor angle Δφ is denoted by ψ_(C). The reference direction is advantageously the direction MOY₁₂ substantially corresponding to the average direction of the bisector between the optical axes O₁Y₁ and O₂Y₂ of the first and second optical sensors D₁, D₂.

In a preferred non-limiting embodiment, the first and the second optical sensors D₁ and D₂ are each formed by a photoelectric receiver such as a photoelectric diode, whose spectral sensitivity can be selected in the spectral range of visible or ultraviolet light or, alternatively, in the range of near or far infrared radiation—for example for nocturnal or all-weather conditions.

The wavelength of maximum sensitivity of the aforementioned photodiodes can thus be chosen in dependence of the specific application of the device for measuring the angular position of a contrasting edge E according to the present invention.

Furthermore, the device comprises a means 1001 (see FIG. 12) for the joint angular vibration of the optical axes O₁Y₁ and O₂Y₂ of the first and second optical sensors D₁ and D₂ in a transverse direction that is different from the given direction.

The given direction, which is substantially parallel to the Z axis, and the transverse direction are different, preferably orthogonal to each other.

Moreover, according to a particularly noteworthy aspect of the device and the method according to the invention, the vibrational law of the optical axes of the first and second optical sensors can be indifferent, aperiodic, unknown and even random without affecting the precision of the angular localization of the contrasting edge.

In FIG. 12, the lens L is stationary and the vibration of the optical axes of the first and second optical sensors D₁ and D₂ is carried out by the translational displacement of a group G formed by the first and second optical sensors D₁ and D₂. Thus, means 1001 for vibrational displacement comprises a support element 1002 which is integral, on the one hand, with said group G formed by the first and second optical sensors D₁ and D₂ and, on the other hand, with a stationary reference mechanical support (not shown in FIG. 12). Means 1001 further comprises a means of translation 1004 which applies to said support element 1002 a stress for controlling the displacement, and hence, for generating a translational displacement in the transverse direction of the group G relative to said stationary reference mechanical support.

Some examples of the means of translation 1004 for generating a translational displacement of the group G are shown and described in patent document FR 2 869 494 (PCT: WO 2005/111536 A1). In these examples, however, the vibrational law is always limited to a periodic vibrational law. Moreover, the device for measuring the angular position of the contrasting edge according to the invention comprises a detection circuit 1020 (see FIG. 12) in order to calculate the angular position ψ_(C) of the contrasting edge E, on the basis of the signals Ph₁ and Ph₂ delivered by the first and the second optical sensors D₁ and D₂, respectively. The detection circuit 1020 is connected to the first and second optical sensors and receives the signals Ph₁(t) and Ph₂(t) delivered by the first and second optical sensors D₁ and D₂, respectively.

The detection circuit 1020 comprises a pair of band-pass filters 1022, each of which is connected to an optical sensor, D₁ or D₂, in order to calculate the differentiated signals, denoted by Ph₁′ and Ph₂′, of the first and second optical sensors, respectively.

Moreover, the detection circuit comprises a calculator 1024 for calculating an output signal Y_(out) and the angular position ψ_(C) of the contrasting edge E depending on the output signal Y_(out).

Furthermore the detection circuit 1020 comprises advantageously a pair of demodulators 1026 connected between the band-pass filters 1022 and the calculator 1024, in order to demodulate the signals delivered by the first and the second optical sensors D₁, D₂ before calculating the output signal Y_(out) and the angular position ψ_(C). The demodulation of these two signals allows one to extract the information of amplitude of these two signals from noise. The amplitude of the differentiated signals is determined using the pair of demodulators 1026. The demodulated signals are then used for the calculation of the output signal Y_(out).

Three techniques of amplitude demodulation, adapted to three different situations, make it possible to extract the amplitude of the two signals in a robust way. A more detailed description of these techniques and the corresponding demodulators 1026 will be given later.

A more detailed description of the method for measuring the angular position of the contrasting edge according to the present invention will now be given.

The lens L placed in front of the two optical sensors confers them a given angular sensitivity. The angular sensitivity function is a bell-shaped (Gaussian-like) function, as shown in FIG. 13. The angular sensitivity functions of the two optical sensors D₁ and D₂ are partially overlapping. The Gaussian-like angular sensitivity function is obtained by slightly defocusing the system consisting of the lens and the optical sensors, whereby the two optical sensors D₁ and D₂ are placed close to the focal plane, between the lens L and its focal plane. The inter-receptor angle Δφ is delimited by the optical axes O₁Y₁ and O₂Y₂ which correspond to the maximum sensitivity of each of the optical sensors D₁ and D₂, respectively.

The first and second optical sensors D₁ and D₂ have each a directional sensitivity, s(ψ) according to the angle Ψ of the light source:

${s(\psi)} = {\frac{1}{\sigma \sqrt{2\pi}}^{- \frac{\psi^{2}}{2\sigma^{2}}}}$

The parameter Δρ is the full width at half maximum (FWHM) of the angular sensitivity function and results directly from the variance σ² of the Gaussian curve by Δρ=σ×2√{square root over (2 in 2)}.

Typical optical characteristics of the optical device are: Δρ=3° and Δφ=2.87°, which gives a ratio Δρ:Δφ of 1.04

Placed in front of a contrasting edge, each optical sensor gives the integral of the product of its angular sensitivity function s(ψ) and the luminance of the contrasting edge E. For example, the contrasting edge is a white/black contrasting edge, where the luminance of the dark zone is 0%, and the luminance of the white zone is 100%. The angular position of the contrasting edge E in the reference frame of the optical sensor is denoted by ψ_(C).

The response of an optical sensor, like a photodiode, is the spatial integral of its angular sensitivity function over the whole white zone:

$\begin{matrix} {{\int_{- x}^{\psi_{C}}{{s(\psi)}\ {\psi}}} = {\int_{- \infty}^{\psi_{C}}{\frac{1}{\sigma \sqrt{2\pi}}^{- \frac{\psi^{2}}{2\sigma^{2}}}\ {\psi}}}} & (1) \end{matrix}$

The integral (Eq. 1) does not have an analytical expression. The error function erf is defined as twice the integral of the Gaussian function with a variance σ²=½ and centred on 0:

$\begin{matrix} {{{erf}(\psi)} = {\frac{2}{\sqrt{\pi}}{\int_{0}^{\psi}{^{- x^{2}}\ {x}}}}} & (2) \end{matrix}$

The primitive of the integral of the angular sensitivity function of an optical sensor is denoted by S, which is equal to 0 in −∞:

$\begin{matrix} {{S\left( \psi_{C} \right)} = {{\int_{- x}^{\psi_{C}}{{s(\psi)}\ {\psi}}} = {\frac{1}{2}\left( {1 + {{erf}\left( \frac{\psi_{C} \times 2\sqrt{\ln \; 2}}{\Delta\rho} \right)}} \right)}}} & (3) \end{matrix}$

Each optical sensor converts the signal into a current. The signals delivered by the first and second optical sensors D₁ and D₂ are denoted by Ph₁ and Ph₂.

Based on Eq. (3), one can calculate the responses of the two optical sensors D₁ and D₂ with their optical axes separated by Δφ and located at

+Δϕ/2 and − Δϕ/2

from MOY₁₂, respectively, as:

$\begin{matrix} {\begin{matrix} {{{Ph}_{1}(\psi)} = {{k \times {S\left( {\psi - \frac{\Delta\phi}{2}} \right)}} + {Offset}_{{Ph}\; 1}}} \\ {= {{k \times \frac{1}{2}\left( {1 + {{erf}\left( \frac{\left( {\psi - \frac{\Delta\phi}{2}} \right) \times 2\sqrt{\ln \; 2}}{\Delta\rho} \right)}} \right)} + {Offset}_{{Ph}\; 1}}} \end{matrix}\begin{matrix} {{{Ph}_{2}(\psi)} = {{k \times {S\left( {\psi + \frac{\Delta\phi}{2}} \right)}} + {Offset}_{{Ph}\; 2}}} \\ {= {{k \times \frac{1}{2}\left( {1 + {{erf}\left( \frac{\left( {\psi + \frac{\Delta\phi}{2}} \right) \times 2\sqrt{\ln \; 2}}{\Delta\rho} \right)}} \right)} + {Offset}_{{Ph}\; 2}}} \end{matrix}} & (4) \end{matrix}$

The gain k represents the difference in luminance between the two surfaces constituting the contrasting edge and depends at the same time on contrast and ambient illumination.

The method according to the invention comprises a step of modulating in a transverse direction, different from the given direction, the amplitude of the signals delivered by the first and the second optical sensors D₁ and D₂, by joint vibration of the optical axes O₁Y₁ and O₂Y₂ of the first and second optical sensors D₁ and D₂, respectively.

The given direction, substantially along the Z axis, and the transverse direction are preferably orthogonal to each other.

The vibration of the optical axes O₁Y₁ and O₂Y₂ of the first and second optical sensors D₁ and D₂, respectively, involves a vibration of the inter-receptor angle Δφ of the first and second optical sensors D₁ and D₂.

Referring to FIG. 11, the vibration of the optical axes O₁Y₁ and O₂Y₂, and therefore of the inter-receptor angle Δφ, is carried out in the plane OXY of FIG. 11.

The step of amplitude modulation results from subjecting the lens L or the group G, formed by the first and second optical sensors D₁ and D₂, to a translational vibration in the transverse direction.

According to a first embodiment of the method and in reference to FIGS. 11 and 12, vibration of the inter-receptor angle Δφ is carried out by relative translational scanning in the transverse direction of the group G formed by the first and second optical sensors D₁ and D₂ with respect to the lens L which is kept stationary with respect to the stationary reference mechanical support. The plane OXY of the FIG. 11 contains this transverse direction, and translational scanning in the transverse direction is denoted by T and is represented by an arrow in FIG. 11.

Thus, a linear displacement of amplitude E applied to the group G formed by the first and the second optical sensors D₁ and D₂, according to the translational scanning T, causes rotation of the inter-receptor angle Δφ and therefore causes each direction of observation OY₁, OY₂ to rotate by an angle Δξ, when the lens L is stationary.

Alternatively, a linear displacement applied to the lens L, of the same amplitudes ε but in the opposite direction, in front of the first and second optical sensors D₁ and D₂ kept stationary with respect to the stationary reference mechanical support, will cause the directions of observation OY₁, OY₂, and therefore the inter-receptor angle Δφ, to rotate by the same angle Δξ.

The temporal evolution of the angular position denoted by Ψ(t) of the contrasting edge relative to the optical axes of the first and second optical sensors D₁ and D₂, as seen from the sensors, is the superposition of the relative angular position of the contrasting edge Ψ_(C)(t) with respect to the bisector MOY₁₂ of the optical axes of the optical sensors, and of the modulation function Ψ_(mod)(t) which results from the vibration effect on the optical axes of the first and second optical sensors:

ψ(t)=ψ_(C)(t)+ψ_(mod)(t)

Each of the signals from the optical sensors D₁ and D₂ is filtered in order to extract their dynamical part by the pair of band-pass filters 1022, which are generally analog circuits (FIG. 11).

Part of the band-pass filters 1022 essentially acts as a temporal differentiator of the signals delivered by the first and second optical sensors D₁ and D₂. The band-pass filters 1022 also cancel the DC component and the high frequency noise, while acting as anti-aliasing filters for the subsequent digital processing. The differentiated signals are denoted by Ph₁′(t) and Ph₂′(t). The expression of the differentiated signals is:

$\begin{matrix} {{{{Ph}_{1}^{\prime}\left( {\psi_{C}(t)} \right)} = {k \times \frac{\psi_{C}^{\prime}(t)}{\sigma \sqrt{2\pi}} \times {\exp \left( {- \frac{\left( {{\psi_{C}(t)} - {{\Delta\phi}/2}} \right)^{2}}{2\sigma^{2}}} \right)}}}{{{Ph}_{2}^{\prime}\left( {\psi_{C}(t)} \right)} = {k \times \frac{\psi_{C}^{\prime}(t)}{\sigma \sqrt{2\pi}} \times {\exp \left( {- \frac{\left( {{\psi_{C}(t)} + {{\Delta\phi}/2}} \right)^{2}}{2\sigma^{2}}} \right)}}}} & (5) \end{matrix}$

Then, the pair of demodulators 1026 isolates the envelope of the differentiated signals.

The output signal Y_(out) is calculated by the calculator 1024 starting from the demodulated signals of Ph₁′(t) and Ph₂′(t). The desired output signal Y_(out) should depend only on the angular position ψ_(C) of the luminance transition zone within the inter-receptor angle Δφ, and should be, in particular, independent of the vibrational law. A convenient expression for the output signal Y_(out) is the difference-to-sum ratio of the amplitudes of these temporally differentiated signals:

${Y_{out}(t)} = \frac{{{{Ph}_{1}^{\prime}(t)}} - {{{Ph}_{2}^{\prime}(t)}}}{{{{Ph}_{1}^{\prime}(t)}} + {{{Ph}_{2}^{\prime}(t)}}}$

This ratio allows for the reduction of the common mode noise originating, e.g., from ambient artificial light (100 Hz and higher harmonics, if the frequency of the mains is 50 Hz), while providing a response which is independent from the luminance, the contrast, and the vibrational law.

Finally, the output signal Y_(out) delivered by the calculator 1024 is:

$\begin{matrix} {{Y_{out}(t)} = {\tanh \left( {\Psi_{c} \cdot \frac{4{{\Delta\phi log}(2)}}{{\Delta\rho}^{2}}} \right)}} & (6) \end{matrix}$

This method to detect a contrasting edge is also applicable to the detection of a bar, a bar being the succession of two contrasting edges of opposite polarities. Indeed, in the same manner, for a bar of width l and centred on ψ_(C), the responses of the two photodiodes are:

$\begin{matrix} {{{{Ph}_{1{bar}}(\psi)} = {{k \times \left\lbrack {{S\left( {\psi - \frac{\Delta\phi}{2} + \frac{l}{2}} \right)} - {S\left( {\psi - \frac{\Delta\phi}{2} - \frac{l}{2}} \right)}} \right\rbrack} + {Offset}_{{Ph}\; 1}}}{{{Ph}_{2{bar}}(\psi)} = {{k \times \left\lbrack {{S\left( {\psi + \frac{\Delta\phi}{2} + \frac{l}{2}} \right)} - {S\left( {\psi + \frac{\Delta\phi}{2} - \frac{l}{2}} \right)}} \right\rbrack} + {Offset}_{{Ph}\; 2}}}} & (7) \end{matrix}$

With this expression, a white bar on a black zone is modelled with a width l>0 and a black bar on a white zone is modelled with a width l<0.

Based on equations 1 and 7, the temporal differentiation realized by the band-pass analog circuit is:

$\begin{matrix} {{{{{Ph}_{1{bar}}^{\prime}(\psi)}(t)} = {k \times \frac{\psi^{\prime}(t)}{\sigma \sqrt{2\pi}} \times {\exp \left( {- \frac{\left( {{\psi (t)} - {{\Delta\phi}/2}} \right)^{2} + {l^{2}/4}}{2\sigma^{2}}} \right)} \times {\sinh \left\lbrack {\left( {{\psi (t)} - \frac{\Delta\phi}{2}} \right) \times \frac{l}{2\sigma^{2}}} \right\rbrack}}}{{{{Ph}_{2{bar}}^{\prime}(\psi)}(t)} = {k \times \frac{\psi^{\prime}(t)}{\sigma \sqrt{2\pi}} \times {\exp \left( {- \frac{\left( {{\psi (t)} + {{\Delta\phi}/2}} \right)^{2} + {l^{2}/4}}{2\sigma^{2}}} \right)} \times {\sinh \left\lbrack {\left( {{\psi (t)} + \frac{\Delta\phi}{2}} \right) \times \frac{l}{2\sigma^{2}}} \right\rbrack}}}} & (8) \end{matrix}$

FIGS. 14 and 15 show the theoretical output signal Y_(out) as a function of the angular position Ψ_(C) of the contrasting pattern, for an edge (FIG. 14) and for a thin bar (FIG. 15). The signal shown in FIG. 14 turns out to be a monotonic and odd function of the angular position Ψ_(C) of the contrasting edge, as expected from the tanh function in Equation. 6. The output signal Y_(out) is symmetrical with respect to the origin (Ψ_(C)=0, Y_(out)(t)=0). In addition, the output signal Y_(out) saturates for the offset positions.

FIG. 15 shows the theoretical response curve to a thin bar. The response curve is again symmetrical with respect to the origin (Ψ_(C)=0, Y_(out)(t)=0) but not completely monotonic since two peaks appear on either side of the abscissa Ψ_(c)=0. It can be shown, however, that these two peaks disappear when the amplitude of the vibration grows larger. In any case, the monotonic and odd central part of the graph can be used for locating a thin bar with high precision. Moreover, it allows one to distinguish between a thin bar and an edge, since in the presence of a bar, the two differentiated signals Ph₁′(t) and Ph₂′(t) have opposite phases, whereas they are in phase in the case of an edge. In addition, as shown in FIG. 16, the ratio of the full width at half maximum (FWHM) Δρ to the inter-receptor angle Δφ can be adjusted such as to make the optical sensor sense and locate only contrasting steps or both contrasting steps and bars.

FIG. 17 shows the experimental signals obtained for the different steps according to the method for optical detection of a rectilinear contrasting edge.

The joint vibration of the optical axes of the optical sensors D₁ and D₂ is carried out by their relative translational displacement (FIG. 12) according to the modulation function Ψ_(mod)(t) represented in FIG. 17 a. Here, a sine wave modulation was chosen for the sake of simplicity. The frequency of this sinusoidal modulation (FIG. 17 a) is 40 Hz and its peak-to-peak amplitude is 0.5°, a value much smaller than the inter-receptor angle Δφ which is 2.87°. FIG. 17 b shows the time course of the actual position of the contrasting edge Ψ_(C)(t) with respect to the bisector MOY₁₂ of the two optical axes O₁Y₁ and O₂Y₂. FIG. 17 c shows the output of the device, which estimates the angular position Ψ(t) of the contrasting edge from both the actual position of the contrasting edge Ψ_(C)(t) (FIG. 17 b) and the modulation function Ψ_(mod)(t) (FIG. 17 a).

FIG. 17 e shows the analogically differentiated signals Ph₁′ and Ph₂′ of the first and second optical sensors D₁ and D₂, respectively.

FIG. 17 f shows the two signals of FIG. 17 e after digital filtering and demodulation.

FIGS. 17 g and 17 h show the difference and the sum of the demodulated signals of FIG. 17 f, respectively.

Finally, FIG. 17 i shows the temporal evolution of the output signal Y_(out)(t), which is the difference-to-sum ratio of the demodulated signals Ph₁′ and Ph₂′. In some parts of the graph (identified in small boxes) it can be seen that the output signal Y_(out)(t) is a linear function of time. Each of these parts is indicative of the presence of a contrasting edge.

FIG. 17 d gives the experimental output signal Y_(out)(t) as a function of the angular position Ψ_(C)(t) of the contrasting edge and compares it with a theoretical curve (dashed line). The theoretical line has been calculated using Equation 6 with Δρ=2.00° and Δφ=2.87°.

In a second embodiment shown in FIG. 18, the means 1001 for vibrational displacement comprises a support element 1006 bearing the lens L. Here, the vibration of the optical axes of the first and second optical sensors D₁ and D₂ is realized by the translational vibration of the support element 1006 and hence of the lens L itself, while the group G formed by the first and second optical sensors D₁ and D₂ remains stationary.

Means 1001 for translational vibration of the lens L comprises a means of translation 1008 for translational displacement of the support element 1006 supporting the lens L. For example, the means of translation 1008 comprises an actuator for generating the translational displacement of the lens L, like the actuator described in patent document FR 2 869 494.

In a third embodiment of the device shown in FIG. 19, the means 1001 for vibrational displacement comprises a support element 1010 integral, on the one hand, with a group G′ formed by the lens L and the first and second optical sensors D₁ and D₂ and, on the other hand, with a stationary reference mechanical support.

Furthermore, means 1001 for vibrational displacement comprises a means of displacement 1012 for applying, to said support element 1010, a stress for controlling the displacement in order to generate a rotational displacement of said group G′ in the transverse direction relative to said stationary reference mechanical support.

In this particular case, according to the method of the invention, a vibration of the optical axes O₁Y₁ and O₂Y₂ of the first and second optical sensors D₁, D₂, in the transverse direction, is carried out by relative rotation of the group G′ formed by the first optical sensor D₁, the second optical sensor D₂ and the lens L, around another axis substantially parallel to the given direction (see FIG. 20).

As an example, referring to FIG. 20, a rotation around the OZ axis by an angle θ applied to the group G′ causes rotation by the same angle θ of the directions of observation OY₁, OY₂ and consequently of the inter-receptor angle Δφ.

Some examples of the means of displacement 1012 for generating a rotational displacement of the group G′ are shown and described in the article by S. Viollet and N. Franceschini in the published article designated A₂: “Visual servo system based on a biologically-inspired scanning sensor”, Sensor fusion and decentralized control in Robotics II, SPIE vol. 3839, pp. 144-155, 1999. In this article, however, the vibrational law is supposed to be exclusively periodic and follows a particular waveform.

A particularly noteworthy aspect of the device according to the invention, is that when the device is implemented as a visual sensor for fixation and/or tracking of a target with a contrasting edge, for example on board of an aerial robot, the natural vibrations caused by the displacement of the robot, or by the turbulences and/or by the robot's engine, will generate vibrations of the platform to which the support element 1010 is firmly attached. The optical axes of the optical sensors D₁ and D₂ will therefore be subjected to the random vibrations of the platform. An advantage of the device according to the invention is its extreme lightness and robustness since additional actuators are not necessary.

A more detailed description of the demodulator 1026, according to the present invention, will now be given with reference to FIGS. 21 to 23.

In the presence of noise, the relevant information contained in the signals Ph₁′ and Ph₂′, that is, the relative amplitude of these two signals, is extracted by demodulation. Several techniques of demodulation exist. Each technique requires adapted filters. In particular, three techniques of amplitude demodulation, adapted to three different situations, make it possible to extract the amplitude of the two signals in a robust way.

In the case considered above, the modulation is sinusoidal (FIG. 17 a) and its frequency f_(a) is a priori known. It is, however, not necessary to know the exact shape of the modulation signal Ψ_(mod). If the modulation signal Ψ_(mod) is not known or not measurable, an asynchronous demodulator can be used. The asynchronous demodulator can also be used when the device for detecting a contrasting edge is mounted on board of a system, such as a machine or a vehicle, of known mechanical resonance frequency. The natural vibration of the system at this frequency will generate signals on the optical sensors which can be extracted with this demodulator. The advantage of this solution is that it does not require any additional sensors such as a rate-gyro sensor to measure the actual vibration, nor any additional actuator to create the scanning process. The realization of the whole device is therefore simple and robust.

FIG. 21 shows an asynchronous demodulator 1026. The front end of the asynchronous demodulator is a pair of selective filters 1030 centred on the frequency f_(a), which extract the differentiated signals, Ph₁′ and Ph₂′, originating from the first and second optical sensors D₁ and D₂, respectively. On each of the two channels, this filter 1030 attenuates all signals whose frequencies are lower and higher than the frequency f_(a), thus improving the signal-to-noise ratio of the signal at frequency f_(a). Furthermore, it comprises a calculator 1032 that generates a sine wave signal and a cosine wave signal at frequency f_(a), which are multiplied by the filtered signals of the two optical sensors. This step results in the translation of the peak from frequency f_(a) respectively to frequency 0 and 2f_(a) in the spectrum. After this step, the continuous component signal is similar to the amplitude of the signal modulated at the frequency f_(a). However, this signal is noisy, because of the frequency 2f_(a) generated by this step.

The asynchronous demodulator comprises moreover means 1034 for calculating the absolute value of the signal in order to relocate the frequency peak 2f_(a) towards 4f_(a). The absolute value is also necessary for the calculation of the output signal Y_(out).

The asynchronous demodulator comprises a low-pass filter 1036 to preserve only the DC part of the signal which is the reflection of the amplitude of the signal with the initial frequency f_(a). This last step of low-pass filtering is more effective to remove the image frequency 2f_(a) due to the frequency doubling resulting from the calculation of the absolute value.

Then, the calculator 1024 (FIG. 12) calculates the difference and the sum of the two demodulated signals, as well as the difference-to-sum ratio which represents the output signal of the device.

According to an alternative and with reference to FIG. 22, the demodulator 1026 is a synchronous demodulator. Such a demodulator is used if the modulation signal Ψ_(mod) is accessible, for example via the control signal of the actuator (e.g., a piezoelectric actuator or a galvanometer) that makes the optical axes of the optical sensors vibrate, or via a sensor (e.g., a rate gyro sensor) that measures the actual micro-vibration of the optical axes. Thus, the modulation signal Ψ_(mod) can be unspecified and measured via a sensor, or preset via an actuator controlled with a predefined command law (Ψ_(mod)).

The synchronous demodulator comprises a selective filter 1040, a multiplier 1042, a means 1044 to calculate the absolute value of the signal and a low-pass filter 1046, all components being comparable with those of the asynchronous demodulator of FIG. 21.

The selective filter 1040 for f_(a) is adapted to retain only frequencies present within the modulation signal Ψ_(mod).

If the modulation is a sine wave signal at frequency f_(a), the selective filter 1040 tuned to the modulation frequency f_(a) will allow the signals of the photodiodes Ph₁′ and Ph₂′ to be set in phase with the modulation signal Ψ_(mod).

The synchronous demodulator comprises means 1042 for calculating the demodulation, and is connected on the one hand to the selective filter 1040 and on the other hand to means 1044 for calculating the absolute value. The means 1042 uses the sine wave signal Ψ_(mod) directly without the need to differentiate or delay the modulation signal because of the intrinsic properties of the sine function.

The calculation carried out during synchronous demodulation is simpler than that carried out during the asynchronous demodulation. It is the most effective demodulation if the modulation signal (Ψ_(mod) or Ψ′_(mod)) is known or can be estimated. The result of demodulation is optimal when the signals coming from the photodiodes Ph₁′ and Ph₂′ are in phase with the signal used for demodulation. This setting in phase can be done by delaying the modulation signal (Ψ_(mod) or Ψ′_(mod)).

Furthermore, the use of a selective filter 1040 tuned to the modulation frequency f_(a) makes it possible to eliminate all the noise present in the spectrum around f_(a). This filter is thus very efficient at the selected frequency but causes an important delay. An alternative is to use several selective filters to reject the identified noise (e.g., 100 Hz, 200 Hz).

According to another embodiment and with reference to FIG. 23, the demodulator 1026 is an envelope detector that detects the envelope of the two signals Ph₁′(t) and Ph₂′(t). Although this type of demodulator gives the least satisfying result in terms of noise rejection, it is the least constraining one of all three demodulators presented here.

The great advantage of the envelope detector over the asynchronous and synchronous demodulators is that it is able to extract the relevant information for the estimation of the angular position of the contrasting edge E without any knowledge of the vibration of the platform supporting the detection device. This characteristic allows for the realization of a particularly simple detection device where the group G′ consisting of the lens and the optical sensors is mounted directly on the platform. The envelope detector will make it possible to use the visual information resulting from the vibration, irrespective of the vibrations of the platform and the origin of this vibration. The realization of an envelope detector is as simple as effective.

In reference to FIG. 23, the envelope detector comprises a comb filter 1050 to remove all the multiple frequencies of 50 Hz, a means 1054 for the calculation of the absolute value and a low-pass filter 1056. The last two elements are both comparable with those of the asynchronous and synchronous demodulators. Means 1054 for the calculation of the absolute value followed by the low-pass filter 1056 extracts the envelope of the signals. The choice of the cut-off frequency of the low-pass filter 1056 is a compromise between an acceptable noise level in the final measurement, and the dynamics of evolution of the measurement. For example, this cut-off frequency may be fixed at 10 Hz.

By way of example, a device for optical detection of a rectilinear contrasting edge according to the present invention, in which the rotational scanning of the inter-receptor angle Δφ of the optical sensors is obtained by vibrations of two photodiodes placed substantially in the focal plane of a lens with a focal distance of f=8.5 mm, gave the following results:

Device for measuring the angular Parameter position according to the invention Δ φ 2.87° Δρ 3°   Range of detectable contrast from 4.5% à 85% Precision of angular localization 0.01°

The precision of angular localization of the contrasting edge is 0.01° which is amazingly greater (287 times greater) than the inter-receptor angle Δφ. The device according to the invention, although being simple, is therefore characterized by an outstanding level of precision.

According to one embodiment, several devices for measuring the angular position can be arranged in a set of devices, in order to increase the total field of view of the detection device until it is rendered panoramic. The set of devices comprises at least one lens, a plurality of N optical sensors (N being an integer and N≧2) placed in the focal plane of the lens and (N−1) detection circuits as described before. Each pair of consecutive optical sensors is connected to a detection circuit to calculate a plurality of output signal Y_(out(N))(t). The output signal delivered by each detection circuit can thus be sampled for subsequent digital processing. The method as described before is applied to each pair of optical sensors to detect a contrasting edge or a bar.

The plurality of optical sensors can be placed in the focal plane of a single lens, the so-called “camera eye” configuration. Alternatively, one or several pairs of consecutive optical sensors can be placed in the focal plane of one of several lenses. Such a set of devices comprises a plurality of adjacent devices, each consisting of a lens and a plurality of optical sensors with different optical axes, the so-called “compound eye” configuration.

In the compound eye configuration, the devices for measuring the angular position are arranged substantially on a spherical or cylindrical surface surrounding a common centre. This common centre constitutes, for the set of devices, the optical centre of an inter-receptor angle Δφ which is the sum over all inter-receptor angles Δφ_(i) of each of the devices taken together.

Furthermore, this set of devices can be generalized for detecting a two-dimensional target comprising at least two contrasting edges, substantially rectilinear and orthogonal.

In the following example the number of optical sensors is limited to three, D₁, D₂ and D₃, in order to facilitate the description. The two pairs of optical sensors, D₁-D₂ and D₂-D₃, are arranged in two directions substantially orthogonal H and V. Each pair of optical sensors is placed in the focal plane of its own lens or all the optical sensors are placed in the focal plane of a single lens. Each lens and each optical sensor are arranged substantially on a spherical or cylindrical surface surrounding at least one common centre. A first centre constitutes the optical centre of an inter-receptor angle Δφ_(H), in the first direction, and a second centre, which can be identical to the first one, constitutes the optical centre of an inter-receptor angle Δφ_(V), in the second direction, substantially orthogonal to the first direction.

In the same manner as before, each pair of optical sensors is connected to a detection circuit to calculate two output signals Y_(out(H))(t) and Y_(out(V))(t), each output signal depending on the position of one contrasting edge of the target.

Finally, in a non-limiting embodiment, the plurality of optical sensors is arranged according to a matrix of optical sensors, each optical sensor constituting a pixel D_(ij). The rectangular matrix of optical sensors made up in this way forms a two-dimensional retina.

A more detailed description of a particularly advantageous embodiment of the device for detecting a contrasting edge according to the present invention will now be given. In general, the direction in which vibration of the optical axes of the first and the second optical sensors D₁, D₂ is carried out can in principle be any direction. In a non-limiting embodiment, however, when this vibration is carried out with the aim of detecting a single contrasting edge E, this direction can advantageously be perpendicular to the given direction along which the contrasting edge E extends.

Such a method may be facilitated by orienting the vibration direction within a plane substantially corresponding to the focal plane of the lens L.

For this purpose and in accordance with a particular embodiment of the detection device according to the present invention, the detection device advantageously comprises a means for orienting an assembly a₁ formed by the support element, the lens L and the stationary reference mechanical support or an assembly a₂ formed by the support element, the group G of the first and second optical sensors D₁ and D₂ and the stationary reference mechanical support. Consequently, the direction of the vibrational displacement of the lens L or of the group G, and thus of the optical axes of the optical sensors, is oriented in the other direction transverse to said given direction of the substantially rectilinear luminance transition zone. Such means for orientation are described in patent document FR 2 869 494 (PCT: WO 2005/111536 A1).

The detection device for measuring the angular position according to the invention can be implemented in a steering aid system.

A steering aid system and its implementation are described in patent document FR 2 869 494. In the steering aid system according to the present invention, the device for detecting the angular position of the contrasting edge E and the detection device for measuring the angular velocity of the luminance transition zone are replaced by the detection devices as described hereinbefore.

The optical device of the device for measuring the angular position, called a position sensor, and for measuring the angular velocity, called a velocity sensor, can be the same. The same holds for the band-pass filters 8 of the detection circuits which compute the derivative signals of the first and second optical sensors.

The methods for measuring the angular position and the angular velocity can be implemented in series or in parallel. The position sensor vibrates. However the vibration amplitude is too small to influence the measurement of the angular velocity.

The angular position sensor and the angular velocity sensor give complementary information for the automatic flight control systems, allowing the staring aid system to fix and fine track a target. 

1. Method for measuring a time delay between two output signals measured by two adjacent optical sensors, the time delay resulting from a movement of an object moving in a given direction, the method comprising: carrying out a spatial filtering of the signals Ph₁ and Ph₂ delivered by a first and a second optical sensor (D₁, D₂), the optical axes of the first and the second optical sensors (D₁, D₂) forming an inter-receptor angle (Δφ) containing the object; calculating a first order derivative and carrying out a temporal low-pass filtering of the signals delivered by the first and second optical sensors (D₁, D₂) in order to obtain signals Ph₁′ (t) and Ph₂′(t); calculating the second order derivative of the signals delivered by the first and second sensors (D₁, D₂) in order to obtain signals Ph₁″(t) and Ph₂″(t); measuring the delay (Δt) between the signals delivered by the first and second optical sensors (D₁, D₂), the delay depending on the angular velocity of the luminance transition zone of the object; characterized in that the step of measuring the delay (Δt) is a feedback loop based on an estimated delay between the temporally filtered signals delivered by the first and second optical sensors, the feedback loop comprising the steps of: estimating a time delay (Δt(t_(n))) at an instant t_(n) between the temporally filtered signals delivered by the first and second optical sensors; calculating the error Ph_(diff)(t_(n)) at the instant t_(n) between the temporally filtered signal delivered by the second optical sensor Ph₂′(t_(n)) at the instant t_(n) and the temporally filtered signal Ph₁′(t_(n)−Δt(t_(n))) delivered by the first optical sensor and delayed by the estimated delay (Δt(t_(n))); and estimating whether Ph_(diff)(t_(n))≠0 i.e. whether a delay error Δt_(error)(t_(n)) exists; if Ph_(diff)(t_(n))≠0, estimating whether the delay error (Δt_(error)(t_(n))) is smaller or bigger than zero on the basis of the error Ph_(diff)(t_(n)) and the second order derivative (Ph₁″(t), Ph₂″(t)) of the signals delivered by the first and second sensors(D₁, D₂) calculating the evolution of the delay (Δt(t)) at the instant t_(n) on the basis of Δt(t_(n)) and Δt_(error)(t_(n)).
 2. Method according to claim 1, characterized in that, the feedback loop on the estimated delay comprises the steps of: filtering the temporally filtered signals of the first optical sensor and/or the second optical sensor with a derivative filter in order to compute the second-order temporal derivative function Ph₁″(t) and/or Ph₂″(t) of the first optical sensor and/or second optical sensor at an instant t_(n); determining the temporal evolution of the estimated delay (Δt_(n)) on the basis of the error Ph_(diff)(t_(n)) and the second-order temporal derivative function Ph₁″(t) or Ph₂″(t) and updating the estimated delay (Δt); these steps being continuously repeated until the estimated delay error (Δt_(error)(t)) between the delayed signal of the first optical sensor and the signal of the second optical sensor is substantially zero.
 3. Method according to claim 1 or 2, characterized in that the update of the estimated delay is increased or decreased with a constant slope according to the sign of the estimated delay error (Δt_(error)(t_(n))) which is the ratio of the error Ph_(diff)(t_(n)) to the second-order temporal derivative function Ph₁″(t_(n)) or Ph₂″(t_(n)).
 4. Method according to claim 1 or 2, characterized in that the update of the estimated delay is increased or decreased by the product of a function α(t).
 5. Method according to claim 1 or 2, characterized in that the error Ph_(diff)(t) between the delayed signal of the first optical sensor and the signal of the second optical sensor or between the delayed signal of the second optical sensor and the signal of the first optical sensor reflects the accuracy of the time delay estimation.
 6. Method according to claim 4, characterized in that the function α(t) is a constant value between 0 and
 1. 7. Method according to claim 4, characterized in that the function α(t) is proportional to the square of the estimated delay (Δt).
 8. Method according to claim 6, characterized in that the function α(t) is modulated by the function |cos(Δt*2πF)|, where Δt is the estimated delay and F is a frequency present on background signals.
 9. Method according to any one of claims 1 to 8, characterized in that the feedback loop on the estimated delay comprises the step of resetting and initializing an integrator by a threshold method applied to the temporally filtered signals delivered by the first and second optical sensors if the error Ph_(diff)(t) exceeds a predetermined value.
 10. Method according to any one of claims 1 to 9, characterized in that it comprises a step of detecting the direction of the movement of the luminance transition zone comprising the steps of calculating a) Ph_(diff)(t) according to Ph _(diff)(t)(A)=Ph ₂′(t−Δt)−Ph ₁′(t) and b) Ph_(diff)(t) according to Ph _(diff)(t)(B)=Ph ₁′(t−Δt)−Ph ₂′(t) comparing the temporal evolution of Ph_(dirr)(t)(A) obtained according to a) to Ph_(diff)(t)(B) obtained according to b) determining on this basis whether the luminance transition zone moves from the first to the second optical sensor or in the opposite direction.
 11. Device for measuring an angular velocity of a luminance transition zone of an object moving in a given direction, the device comprising at least: a first and a second optical sensor (D₁, D₂) having each an optical axis (O₁Y₁, O₂Y₂) that delimit an inter-receptor angle (Δφ), Δφ containing the object; and a lens (L) in order to carry out a spatial filtering of the signals delivered by the first and second optical sensors (D₁, D₂), the lens (L) having an optical center (O), the first and second optical sensors (D₁, D₂) being placed substantially in the image focal plane of this lens (L), the two lines linking the centers (O₁, O₂) of the first and second optical sensors (D₁, D₂) and the optical center (O) of the lens (L) being defined as average directions of observation of the first and second optical sensors (D₁, D₂); and a detection circuit (6) connected to the first and second optical sensors and comprising a band-pass filter (8) in order to carry out a band-pass temporal filtering of the signals delivered by the first and second optical sensors (D₁, D₂), a means 10 for measuring the delay (Δt) between the signals delivered by the first and the second optical sensors, the delay depending on the angular velocity of the light transition zone, and a calculator (12) for calculating the angular velocity starting from the delay (Δt), the angular velocity being the ratio of the inter-receptor angle (Δφ) to the delay (Δt); characterized in that means (10) for measuring the delay (Δt) comprises: a means to sample the temporally filtered signals of the first and second optical sensors, this means following the band-pass filter (8); a means to estimate a delay (Δt(t_(n))) at an instant t_(n) between the temporally filtered signals delivered by the first and second optical sensors, a means to compute the error Ph_(diff)(t_(n)) at the instant t_(n) between the temporally filtered signal Ph′₂(t_(n)) delivered by the second optical sensor at the instant t_(n) and the temporally filtered signal delivered by the first optical sensor and delayed by the estimated time delay Ph₁′(t_(n)−Δt(t_(n))) at the instant t_(n); and a means to compute the evolution of the estimated delay (Δt(t)) starting from the estimated delay (Δt(t_(n))) at the instant t_(n) and the error Ph_(diff)(t_(n)) calculated at the instant t_(n).
 12. Device according to claim 11, characterized in that the first and second optical sensors (D₁, D₂) are each formed by a photoelectric sensor.
 13. Device according to claim 11 or 12, characterized in that means (10) for measuring the delay (Δt) comprises a derivative filter (17) in order to calculate the second-order temporal derivative function Ph₁″(t) and/or Ph₂″(t) of the first optical sensor and/or of the second optical sensor at an instant t_(n) and a means (19) to determine the temporal evolution of the delay Δt on the basis of the error Ph_(diff)(t) and the second-order temporal derivative function Ph₁″(t) and/or Ph₂″(t) and in order to update the estimated delay.
 14. Device according to any one of claims 11 to 13, characterized in that the detection circuit (6) is a bidirectional detection circuit (200) adapted to detect the direction of the motion of the luminance transition zone, that is, from the first to the second optical sensor or from the second to the first optical sensor, the bidirectional detection circuit (200) comprising two means (10A) and (10B) for measuring the delay (Δt) between the signals
 15. Set of devices comprising at least one lens, a plurality of at least three optical sensors and a plurality of detection circuits, characterized in that two adjacent optical sensors form a pair and are placed in the focal plane of one of the at least one lens and connected to one of the plurality of detection circuits, in that each lens of the at least one lens and all optical sensors are arranged substantially on a spherical or cylindrical surface surrounding a common centre, the common centre constituting, the optical centre of an inter-receptor angle which is the sum over all inter-receptor angles of all elementary devices taken together, and in that the devices are devices according to at least one of claims 11 to
 14. 16. Set of devices according to claim 15, characterized in that the at least three optical sensors are arranged in two directions substantially orthogonal.
 17. Set of devices according to claim 15 or 16, characterized in that the optical sensors have a spectral sensitivity which can be the same or different for the individual optical sensors.
 18. A steering aid system for the visual fixation and fine tracking of a target, the target being an object comprising at least one contrasting edge having a luminance transition zone, and for controlling the speed of an aircraft, characterized in that it comprises: a device for detecting the angular position of the contrasting edge; a means to maintain the line of sight of the steering aid system constantly on the position of the contrast edge (E); a device or a set of devices for measuring an angular velocity of the luminance transition zone according to any one of claims 11 to 17; and a means to control the speed of the aircraft according to the angular position and angular velocity of the target. 